Intraprotein electrostatics derived from first principles: Divide-and-conquer approaches for QM/MM calculations

Two divide-and-conquer (DAQ) approaches for building multipole-based molecular electrostatic potentials of proteins are presented and evaluated for use in QM/MM calculations. One approach is a further development of the neutralization method of Bellido and Rullmann (J Comput Chem 1989, 10, 479-487) while the other is based on removing part of the electron density before performing the multipole expansion. Both methods create systems with integer charges without using charge renormalization. To determine their performance in terms of location of cuts and distance to QM region, the new DAQ approaches are tested in calculations of the proton affinity of N(zeta) of Lys55 in the inhibitor turkey ovomucoid third domain. Finally, the two methods are used to build a variety of MM regions, applied to calculations of the pK(a) of Lys55, and compared to other computational methodologies in which force field charges are employed.